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ABSTRACT 

We report a tentative detection of spin correlations in the Tully catalogue of nearby 
galaxies. We define a simple but nontrivial spin correlation function, and find an an- 
alytic estimate of it in the frame of the linear perturbation theory. Then, we present 
the observed spin correlation signal from the Tully galaxies with error bars. The three 
dimensional spin correlation turns out to be significant at the 97% confidence level, 
detected out to a few h~ l Mpc. This observed correlation is consistent with the the- 
oretical prediction based on the gravitational instability picture of galaxy formation. 
An analysis of systematic errors is also presented. The observed strength of correlation 
may be sufficient to significantly affect blank field of weak lensing searches. 

Subject headings: galaxies:statistics — large-scale structure of universe 



1. INTRODUCTION 

Considerable effort has been devoted to the observational detection of any preferred direction 
in the orientation of galaxy spin vectors in the last century. Although several observational stud- 
ies with various sample sizes and methods have reported so far the existence of some pattern of 
galaxy alignments (Strom & Strom 1978; Binggeli 1982 ; Flin 1988; Kashikawa & Okamura 1992), 
there is no established evidence for significant preferential galaxy alignments yet (Dekel 1985; 
Hoffman et al. 1989; Han et al. 1995; Cabanela & Dickey 1999). However, it has been recently 
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pointed out that the sample size of galaxies in the past efforts were too small to detect weak galaxy 
alignments unambiguously (Cabanela & Dickey 1999). 

Theoretically, galaxy spin alignments is also an interesting issue since its existence has been 
predicted by the standard theory addressing the origin of galaxy spins. Very recently, intrinsic 
correlation of galaxy alignments has attracted sharp considerable attention mainly since it could 
contribute as a significant contaminant in blank fields of the statistical search for weak gravitational 
lensing (Catelan et al. 2000; Heavens et al. 2000; Croft & Metzler 2000). In addition, detection of 
the intrinsic galaxy alignments by itself could have a fundamental impact on cosmology since it 
might shed a light on the reconstruction of the initial density field through the galaxy spin field 
(Lee & Pen 2000a, hereafter LPOOa). 

This Letter is aimed for two goals: to estimate the strength of intrinsic galaxy spin correlations 
from a physical theory and to measure the correlation signal directly from observed galaxies. 

2. PHYSICAL ANALYSIS 

In the gravitational instability picture of galaxy formation, the angular momentum of a galaxy 
originates from the local shear due to the proto-galactic tidal interaction with the surrounding mat- 
ter (Hoyle 1949; Peebles 1969). Using the lst-order linear theory (White 1984; Catelan & Theuns 1996) 
(its prediction on the galaxy spin axis has been shown to be quite valid even in nonlinear regime by 
LPOOa), LPOOa have found a general quadratic expression for the expected unit galaxy spins given 
the intrinsic local shears: (LjLj|T) = (1 + a)5ij/3 — aTikTkj where L is a unit galaxy spin, T is a 
unit traceless shear tensor, and a is a correlation parameter introduced by LPOOa to measure the 
strength of the correlation between L and T. If galaxy spins are random, a is zero. While if spins 
are maximally correlated with the shears (Catelan & Theuns 1996), a is 3/5 (it was mistakenly 
cited as unity in LPOOa, see Appendix A in Lee & Pen 2000b, hereafter LPOOb). The real value of 
a must be determined empirically. LPOOa found a = 0.24 in their numerical simulations. 

The above expression of LiLj given T with nonzero value of a implies the existence of the 
spin-spin correlation due to the spatial correlation of the local shears. Here we define a 2nd 
order spin-spin correlation function by rf(r) = (|L(x) • L(x + r)| 2 ) — r/o where r/o is the value of 
(|L(xi) • L(x2)| 2 ) for the case of random spins. r/o = 1/3 for the three dimensional (3D) spins while 
770 = 1/2 for the 2D. Note that |L(xi) • L(x2)| 2 is just cos 2 (A#) where A6 is the angle between the 
two unit spin vectors, and 770 is the expectation value of cos 2 (#) for a random angle 9. Using only 
directional information, the correlation in angle is most readily measured this way. Since we do 
not measure the chirality of each galaxy's rotation axis, the statistic must be symmetric under 180 
degree rotation. 

Using LiLj = (1 + a)5ij/3 — aTikTkj, on e can estimate (|L(xi) • L(x2)| 2 ) = {LiL^LjL'-) such 
that (LiL'tLjfy = (LiLj LiLj} = ((1 + a)<%/3 - af ik f kj }((l + a)<%/3 - af'^} = 1/3 - a 2 /3 + 
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a 2 (TikTkjT^flj) . To calculate {T^T^jT^f'^) , we approximate the unit traceless shear tensor, T as 
a Gaussian variable of the traceless shear, and apply the Wick theorem to derive (T^T^jT^Tlj) ~ 
{f lk f kj ){f' a fl 3 ) + (f^)(f fej f^.) + {f lk fl 3 ){f kj f' a ) = 1/3 + e R (r)/6- Here £ R (r) is the correlation 
of the top-hat convolved density field normalized by £r(0) = 1. 

Hence, we finally have 

r?(r) « ACl(r) (1) 

with A = a 2 /6 for the 3D case. Note that this theoretical prediction has no free parameter and 
follows directly from first principles and an N-body normalization of the shear-spin correlation 
amplitude of a. The approximation is made only when we apply the Wick Theorem. In fact, we 
have verified in numerical realization of a Gaussian random field that the Wick expansion on T 
approximates the true correlation within a few percent relative error. For the case of 2D, it can be 
shown with the same manner that A = 25a 2 /96 . For a detailed derivation of eq. [1], see Appendix 
H in LPOOb. 

On galaxy scales, the power spectrum is well approximated by a power-law, P(k) = k n with 
n = —2, for which £r(v) ~ r _1 . Thus, eq.[l] predicts that for neighboring galaxies (a few h~ l Mpc), 
the spin-spin correlation is of order 1%, decreasing rapidly as r~ 2 . It will decrease more rapidly 
as one moves to larger scales. Thus, we expect negligible spin correlation signal at separations of 
more than a few h~ l Mpc. 



3. OBSERVED SIGNAL 

For this study, we use a subsample of 12122 spiral galaxies from the Tully catalogue (Tully 
2000). The spiral galaxies are selected here as galaxies having a type 0-9 as listed in Third Reference 
Catalogue of Bright Galaxies (de Vaucouleurs et al. 1991, RC3). It is a all-sky low-redshift sample 
with median redshift of 5180 km/s. Among data, we retain necessary five pieces of information for 
each spiral galaxy: the position angle (a), axial ratio (TZ), the equatorial coordinates (1950 equinox 
for RA and DEC), radial distance (redshift), and absolute diameter. The information of the position 
angle and axial ratio comes from the Uppsala General Catalogue (UGC) in the northern hemisphere 
(Nilson 1974) and the ESO-UGC Catalogue in the southern hemisphere (Lauberts 1982). 

In a local spherical coordinate frame (the radial direction is given as the line of sight of a 
galaxy), the spin vector of a spiral galaxy can be determined through a and 1Z by \L r \ = TZ, 
Lg/L^ = tan(-7r — a) upto the fourfold degeneracy: L r — > — L r and {Lg,L^} — > — {Lg^L^}. Since 
we deal with the spatial orientation of the spin axis, only the twofold degeneracy in L r is relevant 
to us. The spin vector expressed in this local spherical coordinates can be easily transformed into 
the global Cartesian coordinates using the given galaxy positions. 

But before determining the spin axis of each spiral galaxy, we consider several possible sources 
of systematic errors which might be involved in measurements of the given a and 1Z. For example, 
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deviation of the shape of spiral galaxies from perfect ellipses, finite thickness of galaxies, differences 
between the UGC and the ESO-UGC catalogue, and others. 

We redistribute the given a and TZ uniformly according to their original order (TZ is redis- 
tributed within each galaxy Hubble type). We then have a perfectly uniformly distributed in the 
interval [0, tt) and TZ in [0,1). The redistribution of TZ statistically removes finite disk thickness 
effects if each disk has the same thickness within each Hubble type. This uniform redistribution 
reduces the final spin-spin correlation strength very slightly, indicating that we do not expect sys- 
tematic effects of global biases to affect our results. However, we would expect the shape-shape 
correlation of galaxies to result in potential problems (7£-related systematic errors). So we choose 
to measure a statistic which does not depend on such correlations. 

One practical way to calculate the true correlation signal free of ^-related systematic errors is 
to use the projected spins: Instead of using the given values of TZ, we can set TZ = for all galaxies, 
which amounts to projecting the 3D spin vectors onto the plane of the sky (the plane normal to the 
line of sight). Using these projected 2D spins, one can construct a real (but 2D) spin correlation 
function, r\2D as a function of the 3D galaxy separations as follows: Let us consider two galaxies 
with unit spins Lj and Lj at Xj and Xj respectively with a separation r = Xj — Xj, r = r/|r|. The 
2D projected unit spin S for each galaxy can be expressed as Sj = Sj/|Sj|, Sj = Lj — (Lj • Xj)xj, 
and the 2D projected unit separation vector is, yj = yi/|y»| , y« = r — (r • Xj)xj. The 2D spin 
correlations, t]id now can be expressed as the difference in angle, /3, between S and y such that 
V2D = (cos 2 (/?j — (3j)) — 1/2 with cos/3j = Sj • yj. This procedure reduces to |Sj • Sj| 2 — 1/2 in 
the small angle (flat sky) limit, but is also valid for all-sky data used here. Note that the two-fold 
degeneracy, (L r — ► —L r ) is irrelevant for the 2D case where L r = 0. 

To find a 3D spin correlation signal free of the systematic errors involved in TZ, one can consider 
the correlation between the 3D and 2D unit spin vectors such that r/(r) = (|L-S| 2 ) — 1/3 where 
S and L are calculated from the uniformly redistributed a and TZ. Here S is free of the 7£-related 
systematic errors by setting TZ = 0, but L still involves "/^-related systematic errors. Thus, some 
fraction of rj(r) must be the false signal. As a device from which the false signal is to be filtered 
out, one can use the following procedure: Generate random 12122 a's. Using these random q's, 
construct 2D random spins, S'. Then calculate the correlation of r/(r) = (|L • S'| 2 ) — 1/3. Repeat 
this process 1,000 times with different sets of 12122 random a's, and take the average, fj'(r) out of 
the 1000 sets of the correlations. Then fj'(r) represents purely the false signal only, otherwise fj'(r) 
would be almost zero. The effective 3D correlation free of the 7£-related systematic errors is now 
obtained as r]E(r) = rj(r) — ff{r). 

One thing one has to keep in mind for the 3D case, however, is the twofold degeneracy. As 
explained in §2, one cannot determine the sign of the radial component of a galaxy spin. Let us 
assume that we have a galaxy pair whose 3D unit spins and projected 2D unit spins are Lj, Sj and 
Lj, Sj respectively. Here each 3D unit spin is determined upto twofold ambiguity, say {L", L|}, 
{Lj, L b j} (L a and L 6 differ by sign of their radial components). Thus, there are four different ways 
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to correlate the spins of this galaxy pair such that = |L? • Sj| 2 , = |L| • Sj| 2 , = |L" • Sj| 2 , 
/j'j = |Lj • Sj| 2 . The spin correlation of this galaxy pair given the twofold degeneracy is now taken 
as the average of the four combinations: (/?• + + + /jj)/4 — 1/3. We apply this to every 
galaxy pair in calculating rj(r). 

Fig.l plots the resulting r]2D and t]e{v) (the square dots) with error bars, and compares it 
with the theoretical predictions (eq. [1]). The error bars are from the numerical formula which we 
have obtained by the Monte-Carlo method. We have generated iV& samples of 10 6 random spins, 
and calculate the standard deviations of the correlations of the random spins We have found that 
the standard deviations are excellently fitted by the following formula: 0.234/ y/Nb for the 3D, and 
l/\/8Nb for the 2D. For the Tully galaxies, Nb corresponds to the number of galaxy pairs in each 
bin. 

Note that the observed correlation is measured in redshift space, which causes the observed 
signal to look flatter than the intrinsic Lagrangian correlation function defined in real space (the 
dashed line). In order to make a better comparison, we convolve the Lagrangian correlation by a 
ID Gaussian filter with a peculiar- velocity dispersion of a v = 150 km/s to obtain the theoretical 
curve in redshift space (the solid line). Since we restrict ourselves to spiral galaxies, we use a 
value of a v close to the observed dispersion of 126 ± 10 km/sec (Davis et al. 1997) for this class of 
galaxies. We note that this very cold distribution of field spirals implies that galaxies did not move 
far from their initial positions, leaving nearest neighbor relations relatively unchanged from their 
Lagrangian values. 

The ID Gaussian filter instead of a 3D one is used to account for the non-uniform galaxy 
distribution. Since galaxies are not uniformly distributed, one has to weight rj(r) by the galaxy- 
galaxy correlation which is proportional r~ L8 when convolving it. This effectively cancels out the 
3D volume factor of r 2 . For the theoretical curves, a power-law spectrum of P(k) = k~ 2 and the 
correlation parameter of a = 0.24 are used. The spin-spin correlation is of much smaller amplitude 
of a 2 £(?") 2 /£(0) 2 than a, and thus is much harder to measure in simulations. Croft & Metzler (2000) 
and Heavens et al. (2000) have attempted it and found signals of similar amplitudes. But due to 
their different statistics from ours, it is difficult to convert their findings into ours. 

The observed correlation signal is in fairly good accord with the theoretical predictions. The 
total significances of the observed 2D and the 3D correlations at the first three bins are at detectable 
95% and 97% confidence levels respectively within a distance of a few /i _1 Mpc, while the correlations 
diminish to zero at larger separations as the theory predicts (see §2). For the total significance, we 
evaluate the \ 2 from the first three bins (degrees of freedom, no = 3) with the null hypothesis of 
random spins. Then, we calculate the probability of x 2 exceeding the observed value (x%d = 7.78, 
x!d = 8.95) and show that the null hypothesis of random spins, i.e., no correlation hypothesis 
at the first three bins is rejected at 95% and 97% confidence levels for the 2D and the 3D cases 
respectively. For this x 2 distribution, the 2D and the 3D correlations are significant at 2.0a and 
2.4a level respectively. 
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Fig. 1. — The spin-spin correlation as a function of galaxy separations, r = cz (km/s). The power- 
law spectrum, P(k) = k~ 2 and the correlation parameter, a = 0.24 are used for the theoretical 
curves. The solid theoretical curve is a Gaussian convolved version of the dashed one with a 
peculiar-velocity dispersion of a v = 150 km/s. 
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Another concerns that could arise in measuring the spin correlation signal is the overlap effect. 
When a separation of a pair galaxy is less than the physical size of the galaxy pair, then the 
spin correlation signal calculated from that close pair might be affected by the overlap effect. We 
remeasure the signal excluding those close pairs whose separations are less than twice the sum of 
the individual radii of the galaxy pair (total 279 galaxy pairs excluded) . But it turns out that the 
excision to twice the radius at which overlap is expected has a negligible effect on the signal (the 
total significance of the observed signal is slightly lowered from 97% to 95% for the 3D, from 95% 
to 94% for the 2D). 

It is also interesting to measure the covariance of rj(r), cov[r](ri),rj(rj)] = (v( r i)il( r j)) /(°~r)i a Vj) 
between the bins using the 1,000 sample sets of random spins. It provides us the degree of indepen- 
dence of the statistics at each bin. We have found cov 

iv(ri),v(rj)] ~ O(10~ 3 ) < 1 for both the 2D 
and 3D cases, suggesting the statistics at each bin is quite independent of that of the neighboring 
bins. 

4. DISCUSSIONS AND CONCLUSIONS 

There are other sources of shear which could change the apparent direction of the spin axis. 
Above all, shears from weak lensing can be a good example. But the shear effect on the galaxy-spin 
axis from the weak lensing is much smaller in this low redshift sample: To estimate the magnitude 
of the weak lensing shear, one has to consider the spatial correlation of the convergence, k defined 
as the ratio of the observed galaxy surface density to the critical surface density (the surface density 
which can just refocus the light beam in the case of null shear). But the expected weak lensing shear 
is less than a percent at 1° for sources at redshift z ~ 1 and grows as n oc z os (Jain h Seljak 1997), 
so the expected weak lensing shear effect on the galaxy spin axis for our sample, z « 0.02 is 0(1O~ 3 ) 
at a separation of 1 ft, _1 Mpc, one order of magnitude smaller than that of the local cosmic shear 
we measure. 

Given that the dominant shear effect on the galaxy alignments at low redshift is from the 
local cosmic shear, the detection of galaxy alignments could affect recent observational surveys for 
the weak lensing effect on the distortion of the galaxy shapes (Van Waerbeke et al. 2000) and the 
interpretation of the lensing effect (Catelan et al. 2000; Heavens et al. 2000). It is beyond the scope 
of this letter to derive a detailed quantitative calculation of the impact of the observed shape-shape 
correlation on weak lensing surveys, which is being addressed in a separate paper (Crittenden et al. 
2000, in preparation). Here we will make a brief qualitative estimate which shows that the intrinsic 
alignments may introduce a non negligible bias. 

To estimate the impact on weak lensing, we project the 3D spin-spin correlation (eq. [1]) 
using the small angle approximation onto the sky to obtain a projected 2D angle-angle correlation 



-8- 



weighted by the galaxy-galaxy correlation function [£(r) ~ (r/Vo) 1-8 with tq = 5.5/i 1 Mpc] : 

V2{y) = t • (2) 

L is the depth of the galaxy distribution in the source plane. The dilution effect causes 772 oc 
1/L, and the usual approximation is to make L — ► 00, stating that two galaxies that appear 
close in projection are not actually close in real space. The fact that galaxies cluster to each 
other modifies this assumption, such that galaxies that are close in projection actually have an 
intrinsic probability of being neighbors that is enhanced by the correlation function. To estimate 
the effect, we use L =1 Gpc (comoving) and model r]2D = 0.01/r for r > 1 Mpc and r\2D = 0.01 
at smaller separations. Placing the source at a typical comoving angular diameter distance of 
2800/i~ 1 Mpc (corresponding to z ~ 1 in a = -65, h = 0.7 universe, see Pen 1999) and assuming 
non-evolving comoving correlation functions, we find an apparent correlation at one arcminute of 
772(1') ~ 0.00024, dropping as 1/y at larger scales. 

Weak lensing surveys measure the correlation function of the shear 73 at separation r, where 
7a, % = 1,2, are the two components of the symmetric traceless shear tensor. Since equation 
(2) correlates galaxies assuming unit ellipticities, it gives a typical intrinsic correlation induced 
shear amplitude c(r) = (71 (0)71 (r) + 72(0)72^)) = (e(O)e(r) cos 2A0) = 2(e) 2 (cos 2 A0 - 1/2) = 
2(e) 2 772D(0 ~ 2 x 10~ 4 , where we used (e 2 ) 1//2 = 0.6 as given by Kaiser (1998). We note that 
recent detection on one arcminute scale (van Waerbeke et al. 2000) finds a candidate weak lensing 
signal for a Top-Hat convolved shear variance, which gives c(l') ~ 5 x 10 -4 when translated 
into our correlation function. This is somewhat larger, but still of comparable magnitude and 
angular dependence as our intrinsic alignment estimate. This model has made many simplifying 
assumptions, most crucial one of which is that the shape-shape and galaxy-galaxy correlation 
functions do not evolve, and that the depth distribution of source galaxies is 1 Gpc, all of which 
may be uncertain by at least a factor of 2. A direct measurement of each of these quantities is 
essential to fully quantify the impact of shape-shape alignment on weak lensing. We note, however, 
that one test can be made which would probe for the presence of intrinsic correlations. The 
observed shear field 7^ can be decomposed into an "electric" (E), or potential flow component, 
and a "magnetic" (B) or curl component (Pen 2000). Weak lensing predicts B = 0, while this is 
not in general the case for shape-shape alignment, which predicts a comparable strength of E and 
B signal. A measurement of B would thus provide a general control experiment. 

We have estimated the expected strength of the spin correlation normalized to numerical 
simulation results (LPOOa) in the frame of the 1st order linear perturbation theory. Then we 
have measured the spin correlation signal directly from the observed spiral galaxies from the Tully 
catalogue, and tested our theoretical predictions against the observed signal. Both the predicted 
and observed amplitudes are of order of 1% at a separation of 1 fo _1 Mpc, in agreement with each 
other. Given our results, we conclude that past surveys would not have had enough sensitivity to 
detect this signal. 
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Qualitative estimates show that these correlations may have a nonnegligible effect on statistical 
weak lensing surveys. We have proposed using the "magnetic" , or curl, component of the apparent 
shape-shape alignment to discriminate between intrinsic shear and lensing induced shear. 
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